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Abstract 

The concept of centrality in high-energy nuclear collisions has recently become a subject of an active 
debate. In particular, the experimental methods to determine the centrality that have given reasonable 
results for many observables in high-energy lead-lead collisions at the LHC have led to surprising behaviour 
in the case of proton-lead collisions. In this letter, we discuss the possibility to calibrate the experimental 
determination of centrality by asymmetries caused by mutually different spatial distributions of protons 
and neutrons inside the nuclei — a well-known phenomenon in nuclear physics known as the neutron-skin 
effect. 


1 Introduction 

In high-energy heavy-ion experiments mm like those now performed at the Large Hadron Collider (LHC), 
the collisions are often categorized according to their centrality aiming to separate the central head-on 
collisions from the peripheral ones in which only the edges of the nuclei collide. This has been generally 
realized by sorting the events according to the amount of particles or energy deposited in specific parts of 
the detector, details varying from one experiment to another laiaEiE]. In its simplicity, the idea is that 
in increasingly central collisions the colliding nuclei disintegrate more completely thereby producing more 
particles. In nucleus-nucleus collisions, an existence of a correspondence between the intuitive geometri¬ 
cal notion of centrality and its experimental determination is supported by the systematics of azimuthal 
anisotropies [TJ [8] in the spectra of low-transverse-momentum (low-p-p) particles [6l [H flOl flT] which, in 
models based on classical collective dynamics [la [la [H] are readily interpreted as reflecting the initial 
geometry of the collision zone. 

Similar experimental procedures in proton-nucleus collisions at the LHC have revealed much stronger 
centrality dependence of hard-process observables, like high-pr jets [T5l [16], than what was anticipated 
e.g. from models for impact-parameter dependent nuclear effects in parton distribution functions (PDFs) 
nattH]. At the same time, however, the minimum-bias versions of the same jet observables — obtained by 
“integrating out” the variable used for the centrality classification — are in good agreement usKHKini with 
the predictions based on collinearly factorized Quantum Chromo Dynamics (QCD). This appears to indicate 
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that the current experimental methods to realize the centrality classification seriously interfere (23 El] with 
the hard processes and, among other proposals [221 1231 (Ml 125]) h has been argued [2S] that even the 
standard energy-momentum conservation plays a significant role. On top of this, the presence of event-by- 
event fluctuations in the initial profile of nucleons inside nucleus can further distort the multiplicity-based 
experimental centrality classification. In the aggregate, the way that the centrality-selected measurements 
in proton-nucleus collisions should be interpreted has turned out largely ambiguous. In nucleus-nucleus 
collisions interferences between the hard processes and the centrality categorization are of less importance 
as the multiplicity of low-pT’ particles used in the events’ centrality classification is much larger and the 
correlations get diluted. 

In this letter, we will discuss a centrality-dependent effect (in its geometric meaning) which stems from 
the fact that in spherical, neutron-rich nuclei the concentration of neutrons is known to increase towards 
the nuclear surface (271 EEj. We demonstrate how this causes effects in electroweak processes that should 
be large enough to be measured at the LHC and could thereby help in resolving open issues concerning the 
relationship between theoretical concept of centrality and its experimental counterpart. 


2 Collision Geometry and Glauber Modeling 


The density of nucleons i in a spherical nucleus A is often parametrized using the two-parameter Fermi 
(2pF) distribution as 





pr/(l + e~), 


( 1 ) 


where the half-density parameter d* controls size of the nuclear core and a* the thickness of the nuclear 
surface. The saturation densities Pq^ are determined by requiring the total amount of protons and neutrons 
to remain constant (we consider only ^^®Pb nucleus in this paper), 


J d^rpP’^(r) = Z = 82, 


j d?Tp^'^{r) = N = 126. 


(2) 


It is an experimental fact [271 EE] that in neutron-rich spherical nuclei the relative amount of neutrons in 
comparison to protons increases near the surface of the nucleus. This is usually referred to as neutron-skin 
effect: In short, the Coulomb barrier that builds up from the positively charged protons limits the extent 
that the proton density can stretch out whereas, being blind to the Coulomb interaction, the neutrons 
can be found further away [29]. The most recent measurement of this phenomenon m for ^^®Pb nucleus 
indicates that the neutron skin does not, unlike its name suggests, have a sharp edge but rather a “halo-like” 
character. As shown in the left-hand panel of Figured] the proton-to-neutron ratio, does not 

drop abruptly but the fall-off towards the edge of the nucleus is gradual. In this plot (and throughout 
the rest of this paper) we have used the parameter values Up = 0.447 fm, dp = 6.680 fm for protons, and 
Un = 0.55 ± 0.03 fm, = 6.70 ± 0.03 fm for neutrons, taken from Ref. |27j . and the error band results from 
adding the variations caused by the quoted two uncertainties in quadrature. 

The Optical Glauber Model m is a commonly used tool in heavy-ion collisions [3I11321133]. In this 
approach, the total inelastic cross section (7“g(s) in collisions of two nuclei A and B with certain center-of- 
mass energy ^/s is given by the integral 




1 _ 
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Figure 1: Left-hand panel: The measured ratio of proton and neutron densities in 208 p]^ g function 
of nuclear radius. Right-hand panel: The ratio Z^siCk)/Nes{Ck) for different centrality classes in p-|-Pb 
(filled yellow rectangles) and Pb-|-Pb collisions (open green rectangles). The heights of the rectangles are 
determined by the uncertainties given for the neutron density in Ref. m- 


where b is the vector between the centers of the colliding nuclei in transverse plane (see e.g. Fig. 20 in 
Ref. [E]), (T™®^(s) is the inelastic nucleon-nucleon cross sectioili], and TabO^) is the nuclear overlap function 


/ OO 

[TP(si) + ri(si)] 

-OO 

[rP(s2) + r|(s2)], 


(4) 


with si ^2 = s ± b/2, and 

/ OO 

dzp^’^ir,z). (5) 

-OO 

The centrality classes Ck are defined as ordered impact-parameter intervals bk < |b| < such that a 

certain percentage {ck+i — Ck)% of the total inelastic cross section accumulates upon integrating, 


1 

(cfc +1 Cfc)% = / d^b 

\ _ g-'rAfl(b)T->' 

^AB -7-00 


6{bk+i-\h\)e{\h\-bk). 


Let us now consider a given hard process (e.g. large-pr direct photon production or equivalent) whose 
contribution to is negligible and does thereby not “interfere” with our centrality categorization. We 
write the cross section for such a process within a given centrality class Ck as 


da%’^iCk) 



^TA(si)r^(s2)dal)-d(Ai?,si,S2), 

id 


(7) 


^The total inelastic cross section should be isospin symmetric as it does not separate the particles of different charge. That 


IS, we assume cr 
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where the sum runs over different combinations of protons and neutrons. The nucleon-nucleon hard-process 
cross sections B, 81 , 82 ) can, in general, depend on the size of the nuclei (via nuclear shadowing or 

equivalent [Ml ESI [36] ) and even on the spatial location of the nucleons inside the nuclei [MllIZI- However, 
to underscore the neutron-vs-proton differences alone we assume here that such effects are approximately 
constant within each centrality class and write the above cross section as 

= {nn)cM“"(Ck) + {ppkMp'iCt) + (»?)&<«<'■*(£») + {pn)c,d4r"(Ck). ( 8 ) 

where 



and refer to nucleon-nucleon cross sections in a given centrality class. It turns out that for 

symmetric A-|-A collisions (for the centrality categories considered here), 

{np)ck = {pn)ck ~ ^ {pp)ct{nn)ck, ( 10 ) 

and we can visualize a given centrality class Ck as simply containing events from collisions of two nuclei with 
effective number of protons Z^{Ck) and neutrons N^{Ck) defined as 

ZesiCk) = ^{pp)ck , = \l{nn)c^ . (11) 

The case of proton-nucleus collisions can be obtained by replacing the 2pF-distribution for the projectile 
proton by a delta function pP = (5(^)(r). As a result, the hard-scattering cross section becomes 

daltiCk) = {p)cJcjl^\Ck) + {n)c,dal^\Ck), ( 12 ) 

with 

rhi~-\-\ 

{i)c, ^ 27r / dbbT^ih), (13) 

and the effective nucleus which the projectile proton “sees” consists of 

ZldiCk) = {p)c, , iVef (Ca:) = (n)c, . (14) 

protons and neutrons, respectively. To evaluate the effective number of nucleons in each case, we have used 
^mei(^ = 2.76TeV) = 65mb (for Pb-I-Pb), and cr“®^(-y/s = 5.02TeV) = 70mb (for p-|-Pb) for the inelastic 
nucleon-nucleon cross sections [38] • The resulting effective proton-to-neutron ratios Z^^^{Ck)/{Ck) and 
Z^Q^^{Ck)/N^Q^^{Ck) computed using Eqs. (fTTI) and (fTH) in various centrality classes are shown in the 
right-hand panel of Figured] While in most central collisions these ratios are very close to the average value 
Z/N = 82/126 ~ 0.65, in very peripheral bins the relative amount of neutrons grows. Since the edges of 
the nuclei are always inside the integration domain in Eq. dOj), the effect of neutron skin in lead-lead starts 
to be visible in more central collision than in the case of proton-lead. This observation also hints that our 
assumption of a point-like proton makes the centrality-dependent effects slightly weaker than what would 
be obtained by assigning the proton with a hnite size. 
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3 Effects of Neutron Skin in Production 


To make the variations in proton-to-neutron ratio visible, an observable for which ^ dapp^'^ ^ 

is required. As the only significant differenc^ between the protons and neutrons is their u and d valence- 
quark content, we need a probe which couples differently to u and d flavours. Here, we will consider 
the production of inclusive charged leptons from —)• decays. This process is theoretically 

particularly well understood, the collinearly factorized perturbative QCD calculations known up to next- 
to-next-to-leading order nulls], and the state-of-the-art calculations incorporate also next-to-leading order 
electroweak corrections [13| on top of this. The existing minimum bias LHC measurements for this process 
in proton-lead and lead-lead collisions are roughly consistent with the pQCD predictions [miiaiie]. 

In the narrow-width approximation, accurate in the asymptotic limit when the decay width Tw of the 
boson is much less than its mass Mw, the leading-order expressions ga HHj for the charged-lepton 
rapidity (y) and transverse momentum (pj-) distribution can be cast as 


d^ 

dydpT 



/ Qiem y 1 

Vsin^^w/ TfwLw 






+ Sq 



X 

/ ± 

< a 

-h 





} 


(15) 


where = 1 ± (1 — and the symbols aem, and Vij denote the fine-structure constant, 

weak-mixing angle and Cabibbo-Kobayashi-Maskawa matrix, respectively. The sum over all partonic flavors 
is restricted by the 6 function which selects only those combinations of quarks and antiquarks for which 
the electric charges e^., eg. sum up correctly to the charge of the lepton. The momentum arguments of the 
PDFs for quarks f^{x,Q^) and anti quarks f^{x,Q‘^) are given by 



M^ey 

2pTy/s 


1 T '^/l ’ 



M^e-y 

2pTy/s 



(16) 


To account for the centrality-dependence of the hard-scattering cross sections, we use PDFs 
for the lead nucleus defined as 


where and are proton and neutron PDFs in a given centrality class Ck, the latter obtained from 
the former based on the isospin symmetry (e.g. = /^’^'^)- Other nuclear effects like shadowing can be 

incorporated as multiplicative correction factors on the proton PDFs, 


In the case of minimum bias collisions these correction factors and their uncertainties have been estimated 
in several global fits gO] [50l [SU [52] . and models for their possible centrality dependence exist [371 El- 
However, in this letter, we will consider only ratios da{i~^)/dcr{i~) and nuclear modifications like this are 
expected to largely cancel even if they were centrality dependent: At a high factorization scale like oc 

^We neglect the photon distribution [ 31 EO] 
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involved here, most of the sea quarks originate from collinear gluon splittings {g —)• qq) which is a flavor- 
independent process for light quarks. Furthermore, the mutually very similar nuclear effects observed in 
charged-lepton [3H [36] and neutrino deep-inelastic scattering [53] on Pb nucleus indicate that the nuclear 
corrections for the valence quarks are also approximately equal, As a result, the 

overall nuclear corrections for and i~ production must be mutually very alike and largely cancel upon 
taking the ratio, separately in each centrality class. In fact, the results shown in plots below have been 
obtained by setting = ff where ff are free proton PDFs for which we have used the general-purpose 
CTIONLO parametrization [54] . Another attractive feature that the ratio dcr{i~^) jda{i~) entails is that many 
experimental systematic uncertainties can be expected to cancel out. Furthermore, one does not need an 
absolute normalization which would need further Glauber modeling. 

The results that follow have been obtained by using the MCFM Monte-Carlo code [551156] at next-to-leading 
order accuracy with the factorization and renormalization scales fixed to Mw. To mimic a realistic 
experimental situation, we integrate over the charged lepton transverse momentum with pT > 25 GeV. 


3.1 Proton-Lead Collisions 


The effect of the neutron skin will be most pronounced in the kinematic region where the large-x nuclear 
valence quarks fEf'‘{x) and are of importance. From Eq. (ITHl) we see that this happens towards 

negative values of y (the “backward” direction). The resulting centrality dependence of ratios da{£~^)/da{i~) 
in p-|-Pb collisions are illustrated in the left-hand panel of Figure [2] by comparing two peripheral classes 70- 
80% and 90-100% to the minimum bias one, 0-100%. Since the main contributions to the cross sections come 
from the ud and du partonic channels m and the sea-quark distributions at small x and large factorization 
scale are approximately flavor independent, f^{x,Q‘^) « f2.{x,Q'^), we can approximate 


^^^P+Pb 

y<o 

fa '\ 


Z!P{Ck)/N^P{Ck) 

fEf\C2,Q^) + f!fy^2^Qy 

^^p^pb 

^a+) 

. 

'z^P{Ck)/N^^\Ck) 



As f^f'^{x,Q'^)/fEf'‘{x,Q‘^) < 1, the derivative of this expression with respect to Z^^^{Ck)/N^Q^{Ck) is 
positive and, in line with Figure [21 the ratio da^^ /da^ decreases towards more peripheral collisions (since 
{Ck) decreases). Towards positive values of y (the “forward” direction), we have 


d^i+Ph 

d(^P+Ph 





( 20 ) 


where we have assumed that the small-x sea quark distributions in Pb nucleus are approximately flavor 
independent at large Q^, Q^) « {x, Q‘^)■ The independence of Z^Q^{Ck)/N^Q^{Ck) explains 

why the centrality dependence of da^^ /da^ virtually disappears towards large values of y. Gurrently, no 
centrality classified proton-lead data for production are available. 


3.2 Lead-Lead Collisions 


The right-hand panel of Figure [2] presents the results in the case of symmetric Pb-|-Pb collisions. As earlier, 
sufficiently far away from the midrapidity, |y| S> 0, we can approximate 


“'^Pb-l-Pb 

da^~ 

“"Pb-I-Pb 


|2/I>0 


U+y [Z™(Cfc)/iVPbPb(c,)] fEf^{x,Q^) + fSf^{x,Q^) 


( 21 ) 
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Figure 2: Left-hand panel: The ratio /da^ in proton-lead collisions for two peripheral centrality class 
normalized by the corresponding ratio in the minimum bias collisions. Right-hand panel: As the left- 
hand panel, but for lead-lead collisions. The small wiggles are residual fluctuations from the Monte-Carlo 
integrations in MCFM. 


where now x = for y <C 0, and x = ior y ':$> 0. Similarly as in the case of proton-lead, the dependence of 

explains the stronger suppression of da{i~^)/da{£~) towards peripheral collisions. In 
comparison to the proton-lead collisions, the effect is better visible even at y = 0 since both the sea-valence 
and valence-sea scatterings depend on while in proton-lead collisions this happens 

only for sea-valence contribution. 

Currently, the most accurate experimental measurements of production in lead-lead collisions come 
from the ATLAS collaboration [l5]. The data for /W~ ratio are plotted in Figure [3] and compared to 
the calculation described in this letter. The data, as well as the MCFM computation, include all events 
within the kinematic region restricted by the lepton transverse momentum pT > 25 GeV, pseudorapidity 
interval < 2.5, missing transverse momentum > 25 GeV, and the transverse mass of the 

system mr > 40 GeV. While our calculation is consistent with the data, x^/^data. ~ 0.6 (V^ata is the 
number of data points), the measurements are clearly not accurate enough and have a too coarse centrality 
categorization to draw conclusions to any direction at this stage. 


4 Summary and Outlook 

We have discussed the generic effects that mutually different spatial distributions of protons and neutrons in 
heavy nuclei are expected to induce in production of W^ bosons at the LHG. The proton density is known 
to fall off more rapidly than that of neutrons towards the surface of neutron-rich nuclei like 208 p|^ which, as 
we have demonstrated, correlates with the sign of the produced W^ boson. Thus, the W^ production could 
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Figure 3: The centrality dependence of W~^IW~ ratio as measured by the ATLAS collaboration |45] com¬ 
pared to the calculation presented in this paper. The reported experimental values of IW~ ratio in each 
centrality class have been normalized to the average one (hF^/hF“)o_8o% = uncertainties have 

been obtained by adding in quadrature the statistical and systematic uncertainties. 


be used to benchmark different experimental definitions of centrality at the LHC. We stress that in this 
paper we have considered the centrality from a purely geometric viewpoint neglecting e.g. possible smearing 
of the effect that could be caused by event-by-event fluctuations. Therefore, our results should be taken as 
hrst, rough estimates of the expected systematics if a given experimental centrality-selection method truly 
reflects the collision geometry. 

The effects caused by neutron skin are, of course, not limited to production but related phenomena 
are to be expected e.g. in the case of high-p-p photons and charged hadrons which can also be measured 
at lower center-of-mass energies like those available at the Relativistic Heavy-Ion Collider (RHIC). In the 
case of lepton-nucleus deeply-inelastic scattering, the neutron skin should affect differently the neutral- and 
charged-current reactions. Thus, the neutron-skin effect could serve as a handle to study the centrality also 
at planned deeply-inelastic scattering experiments like the Electron-Ion Collider |58] or LHeC [59]. 
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